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Abstract 

Fluxon transmission through impurities of different shape in a quasi-one- 
dimensional long Josephson junction is investigated. The junction width is 
significantly less than its length but, at the same time, is of the order of the 
Josephson penetration length Aj or exceeds it. The retrapping current on the 
impurities of the point, hne and rectangular shape is computed as a function 
of the junction width both numerically and analytically. Good agreement 
between the analytic formulae and the numerical simulation results for the 
intermediate (several Aj) junction width is observed. 
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1. Introduction 

The dynamics of magnetic flux propagation in a long Josephson junc- 
tion (LJJ) has been and continues to be a subject of strong theoretical and 
practical interest during the last three decades Q, Isf. The magnetic flux 
quantum in a LJJ is a soliton (also known as fluxon) governed by the well- 
known sine-Gordon (SG) equation. The convenient way to prepare a junction 
with the required properties is to install various inhomogeneities into it. Up 
to now the substantial theoretical work has been devoted to the study of 
the fluxon motion in the one-dimensional (ID) LJJs with point-like 0, H, 0] 
and spatially extended 0, 0, |9| inhomogeneities. Experimental results on 
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the fluxon scattering on impurities are reported in p, Il0| Spatially inho- 



mogeneous Josephson systems with trapped fluxons have been discussed as 



prospective applications, such as fluxon-based information devices [11|, |12 . 

Real LJJs are always two-dimensional (2D), or, more precisely, quasi- 
one-dimensional (QID), in the sense that they have a finite width in the 
direction perpendicular to the direction of fiuxon propagation. Up to now 
the fluxon dynamics in QID junctions has been scarcely investigated as com- 
pared to the pure ID case. Most of attention has been focused on the various 
isotropic 2D structures like oscillons and ring kinks 13|]. It is worth mention- 



ing also the case of window junctions |14l. Il5|. which can be called inverse in 



some sense: it studies point or rectangular junctions embedded in a larger 
two-dimensional superconduncting sample. Quasi-one-dimensional solitons, 
or, in other words, solitonic fronts in infinite (in both x and y directions) 



samples have been studied in detail by Malomed [16|. Several interesting 
results in the absence of dissipation have been reported including the waves 
in the QID sine-Gordon equation, travelling along the soliton line tlTl and 



skyrmion scattering on impurities in the 2D baby Scyrme model [18|. How- 



ever, in the case of the fiuxon dynamics in a LJJ the presence of dissipation 
is unavoidable. 

It is of interest to investigate the QID fiuxon dynamics in the presence 
of spatial inhomogeneities when the junction width is finite. We expect that 
the fiuxon transmission in this case will be significantly enhanced comparing 
to the pure ID case. For example, in the previously studied case of the lattice 
acoustic soliton front interaction with mass impurities [19] it has been shown 
that the front can round the point impurity while a ID lattice soliton gets 
refiected from it. Moreover, the soliton front can overcome even the impurity 
of the infinite mass. To our knowledge the fiuxon interaction with spatial 
inhomogeneities has not been studied except in 16|] , however this paper deals 
only with the infinite sample width. 

Thus, in this Letter we aim at studying the QID fiuxon interaction with 
impurities and finding out how this interaction depends on the junction width 
and other system parameters. In particular, it is planned to find how the 
retrapping current (e.g., the minimal bias current for which the fiuxon prop- 
agation is still possible) depends on the junction width. 

The paper is organized as follows. In the next section, the model is 
described. Section [3] is devoted to the fiuxon transmission through impurities. 
In the last section discussion and conclusions are presented. 
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2. The model 



We consider the QID long Josephson junction (LJJ) subjected to the ex- 
ternal time-independent bias. The main dynamical variable is the difference 
between the phases 6*2 (x, t) — 9i{x, y; t) = y; t) of the macroscopic wave 
functions of the superconducting layers of the junction. The time evolution of 
the phase difference is governed by the perturbed sine-Gordon (SG) equation 
of the form 

(ptt - A0 + [1 + /(x, y)] sin (p = -a(t)t - 7, (1) 



where A0 



'yy 



and the indexes t, x, y stand for the respective partial 



derivatives. In this dimensionless equation the spatial variables x and y are 
normalised to the Josephson penetration depth Aj, the temporal variable t 
is normalised to the inverse Josephson plasma frequency cuj^ [H, H]. The bias 
current 7 is normalised to the critical Josephson current of the junction and 
a is the dimensionless dissipation parameter. The function f[x,y) describes 
the spatial inhomogenity. In the case of point impurities in the general form 
it reads 



N 



fi^^y) = ^l^5ix~ an)5iy - bn) 



(2) 



n=l 



It is supposed that there are impurities in this junction, positioned at the 
points X = an, y = bn, n = 1,2, N, with /i being the "strength" or the 
amplitude of the impurity. Only the microshorts (/i > 0), i.e., a narrow re- 
gions of locally enhanced critical density of the tunnelling superconctucting 
current will be investigated in this article. However, the size an inhomogene- 



ity in experimental samples is finite [10|. Therefore, we consider also the case 



of the line microshort of width dy, stretched along the y direction: 



f{x,y) = ^i6{x) 



y + 



d., 



d„ 



(3) 



Here 6{x) is the Heavyside function. And finally, the rectangular impurity 
of the finite size in both x and y directions 



f{x,y) = II 



"4 



y 



(4) 



will be considered as well. It should be noted that the impurity strength ix 
has different meanings in all three cases (EJ-dl]). For the point impurity ([2]) 
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setting m; = does not automatically yield the pure ID case studied before 
in 0, H, |6[. This case can be retained in the strip impurity case ii dy = w. 
The ID finite-size impurity case 0, [sl, loj is retained in the same way. 

We choose the boundary conditions along the y direction in the von Neu- 
mann form: 

(Py[x,-^,t) = (py(^X,^,t^ =0. (5) 

The boundary conditions along the x axis are periodic: 0(x + L, y; t) = 
(p{x, y; t) + 2-71, where L is the junction length, L > w ^ 1. 



3. Quasi-one-dimensional fluxon transmission through an impurity 

3.1. Numerical simulations 

In order to get an idea about the character of the fiuxon dynamics, the 
numerical integration of the QID SG equation ([T]) has been performed. The 
Josephson phase and its space derivatives are discretized in the following 
way: (l){x,y;t) (j){mh,nh;t) = (pmnit), 

A0 ~ h-"^ (0m+l,n + 0m,n+l + 0m-l,n + 0m,n-l " 4:(f)mn) , whllc the 5-function 

is approximated by the Kronecker (5-symbol. The resulting set of ODE's with 
the boundary conditions were integrated using the 4th order Runge-Kutta 
scheme. Details of the fiuxon propagation through the two identical point 
impurities ([2]) placed along the y axis at Oi = a2 = and bi = —62 = 2 
are given in Fig. [TJ It is important to mention that the dissipation in 
Eq. ([T]) is crucial and the soliton interaction with impurities differs from 
the dissipationless case where the complex resonant behaviour occurs either 



for the 5-like [20| or for the finite-size [21[ obstacles. Far away from the 
impurities the fiuxon exists as only one attractor of the system with the 
velocity, predefined by the damping parameter and external bias. Therefore, 
contrary to the non-dissipative case, the transmission consists of only two 
possible scenarios: trapping and passage. 

For the sake of better visualization the derivative — 0t is plotted on the 
xy plane for the different time moments and for three different dissipation 
values: a = 0.1, a = 0.05 and a = 0.01. Without loss of generality the 
topological charge is assumed to be Q = 1 (soliton) throughout the paper. 
The initial conditions are taken in the form of the approximate (invariant in 
the y direction) soliton solution, placed at the beginning of the junction and 
having with the equilibrium velocity. 
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The fluxon interaction with the rectangular impurity is presented in Fig. 
[2] for a = 0.1 and a = 0.01. In this case, the countour plot of the function 
—(pt is plotted on the xy plane. The fiuxon retrapping current was computed 
numerically as well, showing its steady decay when the junction width was 
increased. These results will be discussed in detail in the next subsection. 

The following conclusions can be drawn from these results. The fiuxon 
shape experiences certain changes after passing the impurities, namely the 
redistribution of the Josephson phase along the fiuxon line in the y direction, 
as well as the slight bending of the fiuxon shape in the same direction. These 
distortions eventually die out after some time, for a = 0.05 and a = 0.1 this 
happens quite soon after the passage through the impurity. For smaller dis- 
sipation {a = 0.01, as shown in Figs. [T}-1) the oscillations of the Josephson 
phase along the y direction seem to survive for much longer time, comparable 
with the fiuxon propagation time along the junction. The numerical simula- 
tions for thinner junctions, w < 10 (not shown in the paper), demonstrate 
that these shape distortions are much wealser and are hardly noticeable. 
Thus, when studying the QID fiuxon interaction with impurities one can as- 
sume with the high degree of certainty that the fiuxon is an almost hard rod 
at least if w ~ Due to finiteness of the junction in the y direction and 

the boundary conditions ([5]), the straight soliton front is the energetically 
most favourable solution and thus it is not possible to observe the arc-lilce 
solitons reported in [l6| . 

3.2. Retrapping current calculation 

Similarly to fiuxon propagation in the ID LJJ 0, Is], |6[, there must exist 
two characteristic values of the bias current, 7c and jthr, Ic > Ithr- Moreover, 
the current- voltage characteristics of the LJJ with impurities have hysteretic 
nature 0, [6|. If 7 > 7c, the pinning on the impurity is not possible and 
there exists only one attractor that corresponds to fiuxon propagation. This 
happens because the bias current is too strong for the fiuxon to get trapped on 
the impurity. In the interval ''jthr < 7 < 7c, at least two attr actors coexist: 
one corresponds to fiuxon pinning on the impurity (there can be several 
different pinned fiuxons if the impurity has finite length ^) and another one 
to fiuxon propagation. If 7 < ''ythri the only possible regime is fiuxon pinning 
on the impurity. The value of 7c is defined only by the properties of the 
impurity, and can be obtained directly from the ID analog. Contrary, for 
the retrapping current the dimensionality of the junction and its width are 
crucial. 
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Far from the impurity, the fluxon kinetic energy is proportional to the 
junction width and equals = 8w[{l — "w^)"^^^ — !]• Here, v^o = [1 + 
(4a/(7r7))^]~^/^ is the equilibrium fluxon velocity in the spatially homo- 
geneous LJJ 0]. In the non-relativistic fluxon case (|foo| <^ 1), one gets 
Voo — 7r7/(4a). By substituting the ansatz 

(l){x,y;t) = 4arctan exp{[2; — X{y, t)] /y/l — v"^} into Eq. ([T]), where X{y,t) 
is the coordinate of the fluxon center of mass, one can obtain the Newtonian 
equation of motion for the fluxon center. 

Since the fluxons under consideration are extended objects in the y di- 
rection, the equation for the center of mass dynamics as well as the impurity 
potential should depend on y. Taking into account the numerical simulation 
in the previous section, we consider the fuxon as an absolutely rigid rod. We 
also mention that the impurity function f{x, y) can be factorised in the cases 
Q-Q, therefore its center of mass dynamics can be effectively projected on 
the X axis and the respective equation of motion can be written as 

mX + maX + = 0, U{X) = -2ti^X + f/o(X) , (6) 

uX 

where the center of mass coordinate X depends only on time. The impurity 
potential Uq{X) now can be calculated from the respective ID problem by 
simply taking away the y-dependent part in Eqs. (l2])-(jl]). The fluxon mass 
m has to be rescaled depending on the type of the impurity. In the pure ID 
case m = 8. This assumption works well only if the impurity is consists of 
lines and/or rectangles. If it has a more complex shape, for example, like 
triangle, the projection on the ID problem becomes more complicated. 

Point impurity. In the point impurity case ([2]), we consider the line of X 
identical equidistant impurities with a„ = 0, n = 1, 2, . . . , which are placed 
symmetrically with respect to the central line |/ = 0. Each impurity creates 
the potential 2/x/ cosh^ X. In order to project the problem on the ID picture 
it is necessary to rescale the fluxon mass from m = 8 to m = 8w/N, thus 
within the kinematic approach the retrapping current can be found as a root 
of the energy balance equation = 2fi, where E^ = m[(l — "w^)"^^^ — 1] ~ 
Awv'^/N + is the fluxon kinetic energy. In the non-relativistic case 

one gets "jthr — {a/TT)^j8^N /w. The correction of the order 0{^) can 
be taken into account with the help of the method, developed in [5|. The 
modification of this method for the QID case is straightforward, therefore we 
describe only the main steps. The improved energy balance relation equates 
the fluxon energy at X = — oo and its losses due to the dissipation, d^E, with 
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the maximal height of the potential barrier U{X): 



IT'-)' 

2^ 



) 



2 



m 



+ /\E = U{Xma.) 



(7) 



where Xmim ^max > Xmin are the extrema of the potential U{X), Xmin — 
— ln(/i/7)/2, Xmax — ^7r7/(2yu). Equation ([7]) is universal and will be used 
for all types of impurities. The energy loss due to the dissipation equals 
AE = 8awN-^ J_^''^{voo - X)dX ~ 41n2y27^a. Inserting this correc- 
tion term AE into the improved energy balance equation and keeping the 
terms up to the order one gets the final corrected expression for the 

retrapping current: 



It appears that the 'jthr{w) relation is obtained simply by dividing the im- 
purity strength by the factor w/N, moreover, the O(a^) correction does not 
depend on the junction width at all. From this expression one can clearly 
see that if only the 0{a) term is taken into account, the retrapping current 
disappears as w — oo, thus, in the infinitely wide junction a fiuxon always 
passes the impurity. The second term in Eq. (|8]) does not depend on w, and, 
therefore it may lead to the wrong conclusion that 'jthr does not tend to zero 
as w — )■ oo. However, it should be noted that this term has been derived 
under the assumption of w being finite. 

The numerical simulations confirm the main suggestion of this Letter: the 
retrapping current decays with the growth of the junction width. The ap- 
proximate expression (|H]) appears to be in good agreement with the numerical 
data for smaller values of w, while for larger w's the numerical and analytic 
data diverge (however, the agreement remains to be satisfactory). All this 
is presented in Fig. [1^ where the analytic results are given by solid lines 
and the numerical data are shown by markers. In the case of two and three 
impurities, placed along the line in the y direction, the retrapping current 
virtually does not depend on the distance ay between them. For instance, 
the markers that correspond to aj, = 2 (x), a^, = 3 (o), = 4 (♦) at = 2 
and w = 5 are almost indistinguishable. 

One may consider the excentrically placed (i.e., lying away from the x = 
axis) impurity. The kinematic approach does not distinguish the impurities, 
placed at the different positions. The numerical simulations show that the 
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difference between the retrapping currents in both the cases is very small. 
For example, 'jthr = 0.03770 for the centrally placed impurity, ai = bi = 0, 
at w = 2 (the rest of the parameters are as in Fig. [1^) and 'jthr = 0.03776 
for the impurity at ai = 0, 6i = w/4 . 

Line impurity. In the case of the impurity line (|3]), the same energy 
balance expression holds as for the point impurity but the fiuxon mass is 
assumed to be m = 8w/dy. The rest of the calculation procedure for the 
retrapping current is the same as for the point impuritiy. As a result, we 
obtain the the final expression for the retrapping current: 



Numerically computed 'fthr appears to be in good correspondence with the 
approximate expression ([9]) as shown in Fig. [3b. Similarly to the point im- 
purity case, the discrepancy between the analytical and numerical results 
increases at larger w. In the limit dy/w — ?■ 1, the effective ID picture is re- 
stored because the impurity strip crosses the whole junction in the y direction. 
Thus, the retrapping current attains the value 'jthr = a(-\/8/I — 4a ln2)/7r for 
the ID soliton case (shown by the thick horizontal line in Fig. |3]d). Thus, we 
observe that if the strip impurity consitutes, for example, about 1/3 of the 
junction width, the retrapping current is about 40% less than the respective 
ID value. 

Rectangular impurity. Finally, we consider the rectangular impurity 
dH). In the point-particle description, the fiuxon "feels" this impurity as 
the potential Uo{X) = 2/i[tanh(X + 4/2) + tanh(4/2 - X)] Q. In this 
case, the energy balance reads 8w[{l — f^)^^^^ ~ 1] = 4/i(ij, tanh((ia;/2), from 
which the expression for the threshold current can be easily computed. For 
large it can be complemented by the correction 0(a;^) computed for the 
ID case in 0, S]. The fiuxon mass should be renormalized as 8 — )■ 8w/dy. 
In this limit it is assumed that the impurity creates the step-like potential 
Uo{X) = 2/x[tanhX -|- 1], which has a minimum at Xmin = — arcsechy^TTT/Ti 
and a maximum at X^ax = arcsech y^vrT/Tl. The energy losses due to the 
dissipation AE = 8awdy^ j^^"'' {voo — X)dX ^ ijLdyW^^a{hi[iJ,/ {tc^)] + 
C(l)} should be substituted in the energy balance equation The final 
approximation reads 





Ithr 



a d, 



71 




4:\/ /itanh(4/2) — a In 



fl w 

dy 



) 



(10) 
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Note that at variance with the point and hne impurities, the second correction 
depends on the both junction and the impurity width. In the ID hmit 
{dy/w — > 1), the formula ffTOl) retains the form, found in 0, Isf. In the limit 
w — )■ oo, the second term decays to zero, but much slower than the decay 
law 0{w~^^'^) of the first term. Once again, the second term is a correct 
approximation only for finite values of w. 

In the opposite limit, — )■ 0, one can use the following approxima- 
tion: tanh(X + dx/2) + tanh((ia; — X) = sinh(i2:/[cosh^ X + sinh."^ {dx /2)] ~ 
dxSech^X + 0{dl). Thus, the effective potential U{X) is approximately the 
same as for the line impurity up to the coefficient sinh dx — dx- Therefore the 
formula ^ can be modified accordingly and finally the retrapping current 
reads 



The results of numerical simulations given in Fig. H] demonstrate satis- 
factory correspondence with the analytical approximations ( !TO|) - (fTT]) . The 
deviations from the numerical results appear to be stronger as compared to 
the point and line impurities, mostly because we were unable to provide the 
0{a'^) correction that works equally well in both limits dx ^ ^ and dx oo 
and also due to higher degree of the fluxon deformation during the interation 
with the impurity. It should be mentioned that even for dx ^ 2) one can 
consider the "long" dx limit and it works fairly well. Indeeed, the markers 
that correspond to dx = 3 (©), dx = 4: (o) and dx = 5 (♦) seem to be almost 
indistinguishable in Fig. HI 

4. Discussion and conclusion 

To summarize, we have shown that quasi-one-dimensional fluxon passage 
across microshorts is significantly enhanced in comparison with the purely 
one-dimensional case. The retrapping threshold current decays with the junc- 
tion width approximately as w~^^'^, according to the kinematic approach. The 
numerical simulations support this dependence for intermediate (several Aj) 
values of w. With the increase of w, the discrepancy between the numerical 
results and the kinematic approximation becomes more distinct. The reason 
of this discrepancy lies in the fact that for w ^ 1 the fluxon cannot be longer 
considered as an completely rigid object and its deformation the y direction 
should be taken into account. 




(11) 
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One can formulate the following simple argument that explains the en- 
hanced fluxon transmission across the obstacle of the width dy in the QID 
case. The obstacle can be described as a localized (both in the x and y di- 
rections) potential barrier. Only the central part (|?/| > dy/2) of the initially 
homogeneous in the y direction fluxon takes part in the interaction process, 
while the marginal areas dy/2 < \y\ < w/2 do not. Thus, if the energy in 
the tails is sufficient enough to overcome the barrier, the fluxon will pass. If 
w — )■ oo, the energy in the non- interacting part of the fluxon tends to inflnity, 
and, consequently, it will overcome any localized obstacle. 

Finally, we would like to outline the future research directions. In our 
opinion they lie beyond the rigid rod approximation used in this Letter. In 
this case the fluxon center of mass depends also on the transverse coordi- 
nate y and the effective equations of motion are PDE's and not ODE's, like 
Eq. ([6]). Accounting for the dependence on y becomes important in the fol- 
lowing cases: (i) the junction width is too large and the flexural oscillations 
along the fluxon line appear; (ii) the dissipation is rather small and/or the 
junction length is L ~ ly, thus the spatial distortions of the fluxon line do 
not have enough time to die out. These cases are also interesting in con 



nection with the arc-like fluxons found in the inflnitely wide junction [16 



and the variety of the excitations that travel along the fluxon crest in the 



dissipationless 2D junctions [17j. The combined effect of the dissipation and 
the finite junction width makes the transverse-invariant fluxon front the en- 
ergetically most favourable solution. However, other attractors of the system 
which are non-transverse-invariant may exist as suggested by Figs. [TJ-1, es- 
pecially when dissipation is small. Finding them and investigating how they 
will manifest themselves on the junction current-voltage dependence is an 
important problem. 
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t=24 



t=64 



t=80 



t=120 




Figure 1: (Colour online) Time evolution of the Josephson phase derivative in the L = 
70 X w = 10 junction with N = 2 impurities with the distance ay — 4 between them at 
H = 0.5 and a = 0.1, 7 = 0.03 [(a)-(d)], a = 0.05, 7 = 0.02 [(e)-(h)], a = 0.01, 7 = 0.0035 
[(i)-(l)]- 
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Figure 2: (Colour online) Contour plot of the time evolution of the Josephson phase 
derivative in the L = 60 x w = 10 junction with the rectangular impurity (|3]) with /i = 0.5, 
3, dj, = 5, a = 0.1, 7 = 0.08 [(a)-(d)] and a = 0.01, 7 = 0.008 [(e)-(h)]. 
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Figure 3: Retrapping current as a function of the junction width w for a = 0.1, /z = 0.5. 
Panel (a) corresponds to the case of point impurities (H)) with sohd lines illustrating the 
approximate formula ^ a.t N = 1 (line 1), = 2 (line 2) and N = 3 (line 3). Markers 
correspond to the results of numerical simulations at iV = 1 (u); N — 2 with ay — 1 (+), 
Cj, = 2 (x), Cj, = 3 (o), Uy — 4 (♦) and ay — 7 (o); = 3 with ay ~ 1 (A) and ay = 2 
(V). Panel (b) corresponds to the case of the line impurity ([3]) with dy — 1 {1 and +) and 
dy = 2 (2 and x). Markers correspond to numerical results and solid lines correspond to 
the formula ([9]) . Thick horizontal line corresponds to the retrapping current on the point 
microshort in the pure ID case. 
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Figure 4: Retrapping current as a function of the junction width w for the rectangular 
impurity (|4|) at a = 0.1, ji = 0.5 and dy = 1. The markers correspond to the cases dx = 1 
(+), dx ~ 2 (x), dx = 3 (®), dx = i (o) and dx = 5 (♦). The analytical approximation 
((TOl) is given by the solid lines 1 — 5 for da; = 1, 2, . . . , 5, respectively. The inset shows the 
retrapping current dependence on dx at dy = 1 (+), dy — 2 (x) and w ~ 3. Two thick 
solid lines correspond to the analytical approximation (jlOp at dj; — > oo for dy — 1 and 
dy = 2, thin solid line corresponds to the approximation at dy = 1. 
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